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ABSTRACT 



Aims. We discuss the spin-orbit orientation of the Fomalhaut planetary system composed of a central A4V star, a debris disk, and a 
recently discovered planetary companion. 

Methods. We use spectrally resolved, near-IR long baseline interferometry to obtain precise spectro-astrometric measurements across 
the Br-y absorption line. The achieved astrometric accuracy of ±3yuas and the spectral resolution R = 1500 from the AMBER/VLTI 
instrument allow us to spatially and spectrally resolve the rotating photosphere. 

Results. We find a position angle PA stal - = 65° ± 3° for the stellar rotation axis, perfectly perpendicular with the literature measurement 
for the disk position angle (PA disk = 156!0 ± 0:3). This is the first time such test can be performed for a debris disk, and in a non- 
eclipsing system. Additionally, our measurements suggest unexpected backward-scattering properties for the circumstellar dust grains. 
Conclusions. Our observations validate the standard scenario for star and planet formation, in which the angular momentum of the 
planetary systems are expected to be collinear with the stellar spins. 

Key words. Stars: individual: Fomalhaut; Stars: planetary systems; Stars: rotation; Methods: observational; Techniques: high angular 
resolution; Techniques: interferometric 
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1. Introduction 

At 7.7 pc from the Sun, Fomalhaut is one of the closest main 
sequence star surrounded by a spatially resolved debris disk. It 
has been studied at many different wavelengths ranging from the 
visible to the sub-millimetre regime. At long wavelengths (far- 
infrared and beyond), the disk shows a pair of intensity maxima 



inter preted as the ansae of an inclined ring a bout 140 AU in ra- 
dius dHolland et alj2003HStapelfeldt et alj2004l) . This particular 
geometry was sub sequently confirmed with HST/ACS imaging 



geometry w as sub sequently confirmed with H5> I /AGs imaging 
dKalas etai]|2005l) . which also provided high accuracy estima- 
tions of the disk inclination and position angle: /disk = 65°9 + 0°4 
and PAdisk = 156°0 + (f.3, under the assumption that the disc 
is intrinsically circular. In the scattered light images, dust parti- 
cles are confined in a narrow ring about 25 AU in width, with 
sharp edges th at suggest the presence of a planetary body at 
about 120 AU (Ouillen 2006). It was also noted that the centre 
of the ring is shifted by 15 AU with respect to the central star 
position, pointing on still poorly understood dynamical effects 
that should also be related to the presence of planetary-mass 
or sub-stellar companions. Recently, the presence of a planetary 
companion at the expected orbital distance (119 AU) was finally 
confirmed with HST/ACS coronagraphic imaging (Kala s et alj 
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120081) . which significantly boosted the general interest in this 
system. 

The circumstellar dust around stars with ages above ~10 Myr 
might be partially a remnant from the primordial disc. The disks 
could also be replenished by the pop ulations of planet esimals 
that were not used to build up planet (Mann et al. 2006). These 
leftovers are supposed to produce dust by mutual collisions or 
cometary activity. The study of those debris disks provides one 
of the best means to explore the properties (size, density, orien- 
tation) and evolution of planetary systems. 

In standard planet formation scenarios, the disk and the as- 
sociated planets are expected to be in the equatorial plane of the 
central star, like our Kuiper Belt which is within 2° of the in- 
variable plane of the Solar system dBrown & Panll2004l) . Yet, as 
far as debris disks (and circumstellar disks in general) are con- 
cerned, this simple observational prediction has never been put 
into test. More generally, the common assumption of collinear- 
ity between planetary orbital axis and stellar spins has only been 
tested on the Solar system, with a relevant accuracy (~ 1°), and 
on a few transiting extrasolar giant planets with a rather poor 
accuracy (~ 10°, by the mean o f the Rossiter-McLaughlin spec- 
troscopic effect, iRossiteij 1 1 9241) . In at least one case, XO-3b, a 
large spin- orbit misalignment is suggested by the observations 
(70° ± 15°. lHebrard et alj|2008l) . Additionally, the recent contro- 
versy about the spin-orbit alignment of HD 17156b illustrates 
well the difficulty of such a techni que, with m e asured spin-orbit 
misalignme nts of 62° ± 25° by iNarita et all d2008) and later 
9°4 + 9°3 bv lCochran et al.1 d2008l) . Up to now, this test has never 
been performed in non-transiting systems, because it is generally 
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Fig. 1. Projected baselines of the observations. Each observation 
(a, b, c and d, draw with different symbols) is composed of 3 
baselines (1,2,3) along which the projections of the 2D spectro- 
astrometric displacements are measured (see Fig. [2]). 



Table 1. Observation log with the number of files (#), the num- 
ber of frames per file (NDIT) and the elementary integration time 
(DIT). Calibration stars are in italic. 



Date 


Baseline 


Target 


# 


\iniT„ hit 


Obs. a) 








lotal: oUUs 


ZOOo-Iz-Oj 100:20 


D0-ri0-(-jl 


Fomalhaut 


J 


zUUxU.zs 


ZUOft- 1 Z-Vj 1 U t . t j 


DO-riO-O 1 


soAqr 


J 


ZUUxU.zs 


2008-1 2-05T01 -30 


DO-HO-G 1 


Pnmalhflnt 


5 


200x0 2s 


2008- 12-05T0 1:48 


D0-H0-G1 


yScl 


5 


150x0.5s 


2008-12-05T02:06 


D0-H0-G1 


Fomalhaut 


5 


200x0. 2s 


2008-12-05T02:24 


D0-H0-G1 


88Aqr 


5 


200x0.2s 


Obs. b) 








Total: 200s 


2008-12-08T02:15 


A0-K0-G1 


Fomalhaut 


5 


200x0.2s 


2008-12-08T02:31 


A0-K0-G1 


88Aqr 


5 


200x0.2s 


Obs. c) 








Total: 200s 


2008-12-09T00:10 


A0-K0-G1 


Fomalhaut 


5 


200x0.2s 


2008-12-09T00:27 


A0-K0-G1 


88Aqr 


5 


200x0.2s 


Obs. d) 








Total: 200s 


2008- 12-09T0 1:31 


A0-K0-G1 


Fomalhaut 


5 


200x0.2s 


2008- 12-09T0 1:51 


A0-K0-G1 


yScl 


5 


120x0.5s 



impossible to recover the orientation of the rotational axis of the 
central star. 

With the advent of spectrally dispersed optical long base- 
line interferometry, it becomes possible to put this hypothesis 
into test by spatially and spectrally resolving the stellar photo- 
spheres, even marginally. Indeed, this technique allows the ori- 
entation of stellar axes to be quickly (geometrically) recovered 
by differential astrometry across a photospheric line enlarged 
by the stellar rotation. The purpose of this Letter is to perform 
this test in the case of Fomalhaut, which advantage ously com- 
bines a fast projected rotational velocity of 93 km/s dRover et al.l 
12007]) and a relativel y large apparent diameter of <£ = 2.2 mas 
dDi Folco et aT]|2004l) . 



2. Observation and data reduction 

Data have been colle cted at the Very Large Telescope 
Interferometer (VLTI, lHag uenauer et al.l 120081) with the 



AMBER beam combiner dPetrov et al. 20071) recording spec- 
trally dispersed fringes between 1.92 - 2.28 yum at a spectral 
resolution of R = 1500. The Auxiliary Telescopes (ATs) were 
placed at stations D0-H0-G1 and A0-K0-G1. Observations 
were assisted by the fringe-tracker FINITO, allowing an indi- 
vidua l exposure time of 200 ms in AMBER dLe Bouquin et alj 
2008). The seeing ranged 0.8" to 1.2". FINITO provided good 
locking ratio and fringe tracking performance. Observations of 
calibration stars were interleaved with those of the science tar- 
get. Unlike most classical interferometric observations, the only 
prerequisite for calibrators is not to show spectro-astrometric 
features around the Br-y line. We therefore favour cool giants 
because they have less pronounced Br-y lines and smaller 
rotation velocities. The observation log is given in Table [TJ 
Regarding the magnitudes, Fomalhaut is 1.00 K mas , y Scl is 
1 .72 K mag , and 88 Aqr is 0.98 K mag . Because it is fainter, data on 
yScl have been obtained with slightly longer integration time. 
This has no incidence on the differential analysis made in this 
paper. 

Phases of interferometric fringes were computed using the 
latest version of the amdlib package (version 2.99, Chelli et 
al., private communication) and the interface provided by the 



Jean-Marie Mariotti CentetQ. Following a standard fringe selec- 
tion criterion, we have kept the 80% best frames. Consecutive 
observations of Fomalhaut were grouped into single data points 
with enhanced signal-to-noise ratio (SNR), referred to as obser- 
vations a, b, c, and d in TableQ] Accordingly, the observations of 
the calibration stars were also appended into single, high quality 
measurements of the instrumental phase. Such an average in- 
troduces baseline smearing. Considering the worst case of 2h, 
we estimated the smearing of the astrometric signal to be less 
than 20%. This is a reasonable price to pay to obtain clear de- 
tections on all dataset. After the calibration of scientific data, we 
therefore end up with a set of 4 observations of Fomalhaut (see 
Table Q]), each of them composed of the 3 phase spectra from 
the AMBER spectrograph. While signal was also detected in the 
visibility, we decided to focus the analysis of this letter on the 
phase, which contains the astrometric quantity. 

The unavoidable atmospheric dispersion appears as a global 
curvature of the fringe phases over the K band. It cancels per- 
fectly when computing the sum of the 3 phase spectra, demon- 
strating that the baselines are consistent to each others in term 
of sign (so in direction on sky). Subtracting the instrumental 
phase from the scientific phase removed all fine spectral features 
down to a level of 0°5 (uncertainty level). However, atmospheric 
dispersion was not perfectly canceled out. It was efficiently re- 
moved by fitting the spectral range 2.03-2.27 /mi with an ad-hoc 
polynomial law of order 3. As a matter of fact, phase measure- 
ments are relative to the continuum value which we suppose to 
have a zero phase. It may not be perfectly true because of the 
gravity darkening in the fast rotating atmosphere. However, such 
effect would not bias our measurement. 

The calibrated differential phases <j> were converted into 
differential astrometric shifts p using the well-kno wn formula 
for m arginally resolved interferometric observation dLachaumd 
120031) : 

d> A 



P 2n B 



(1) 



B is the length of the interferometric baseline and A is the effec- 
tive wavelength of the spectral channel considered, p represents 
the projection, on the baseline direction, of the 2D photocenter 
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Fig. 2. Top: continuum-normalized spectrum obtained with 
AMBER for a spectral window centered on the Br-y 
line. Bottom: Differential phase measurements converted into 
spectro-astrometric shifts p(A) in mas. Colours represent the po- 
sition within the line, from blue to red. The solid lines link the 
results, reprojected on the baselines, of the global 2D astrometric 
solution p(A) performed independently for each spectral chan- 
nel. Phase spectra are shifted vertically by 0.08 mas. 



in the plane of the sky, hereafter called p. A single astrometric 
solution (i.e. a single 2D vector p) was fitted to all observations 
available within a single spectral channel (i.e. the 12 projections 
p). This global fit is perfectly linear and has the advantage of 
showing the ultimate astrometric performance. However, the re- 
sulting astrometric signal p(A) can be slightly biased (towards 
reduced amplitude) because of the smearing that may occur in 
case of different spectral calibration between the observations. 
This has no impact on the astrophysical interpretation made in 
the following. The differential astrometric shifts p(A) are repre- 
sented by the dots and the error bars in Fig. [2] The result of the 
global fit by the 2D vector p(A) is represented, reprojected on the 
baselines, by the solid lines that connect the fitted values. 

Uncertainties were propagated into the astrometric vector p 
by standard formulae. We did not take into account the uncer- 
tainty on the baseline length nor on the spectral calibration be- 
cause they affect all spectral bins in the same way and therefore 
have negligible effect on our final observable. The astrometric 
error ellipses were found to be almost circular thanks to the rela- 
tively uniform u, v plane coverage provided by our observations 



(see Fig.[TJ. Uncertainties are in the range +2 //as to +6 //as, with 
typically +3 //as across the Br-y line. It is among the most pre- 
cise spectro-astrometric measurements ever achieved, if not the 
best. 

When talking about potential bias, our main concern was the 
presence of a possible artifact on the differential phase (p(A) due 
to the particular spectrum of our target (A4V, i.e. strong Br-y 
absorption line) in comparison to the calibrator spectrum (K gi- 
ants, i.e. almost no Br-y line). However, an instrumental artifact 
linked to the AMBER spectrograph would rotate on sky between 
observations a), b), c) and d), following the rotation of the asso- 
ciated baseline triangles (see tableland fig. [TJ. The good agree- 
ment between the signal obtained during 3 different nights and 
on different configurations confirms its astrophysical origin. As 
an additional test, we observed other fast rotating stars and ob- 
tained various kinds of spectro-astrometric signals, as well as 
expected non-detection on small photospheres. 

3. The Fomalhaut rotation axis 

Our observations represent the most precise spectro-astrometric 
displacements ever measured with the AMBER instrument 
(±3 //as in the Br-y line). To extract the scientifically useful in- 
formation, we plot the spectro-astrometric photocenters onto the 
plane of the sky. As a first check, we plot the photocenters in two 
spectral windows of continuum around the Br-y line (Fig. [3^). 
The dispersion of the data points as a function of wavelength in 
the continuum is fully compatible with our estimated uncertain- 
ties (±3 //as) and shows no systematic behaviour or unexpected 
features at our precision level. 

On the contrary, a significant astrometric displacement is de- 
tected when plotting the signal across the Br-y line (Fig. [3p). 
The photocentre is displaced in one direction on the blue wing 
of the line, and is identically displaced on the opposite direction 
on the red wing. The shift is compatible with zero in the centre 
of the line, and its maximum amplitude in the line wings reaches 
about 15//as. This is a clear detection of the effect of rotation 
on the stellar photosphere, from which the Br-y absorption line 
originates. Interpreting such a small astrometric shift does not 
require deep modelling of the fast rotating atmosphere. In par- 
ticular, the astrometric displacement is necessarily perpendicular 
to the position angle of the stellar rotation axis. Fitting our com- 
plete data set, we find the stellar rotation axis to have position an- 
gle PA stal = 65° + 3°. This value is perfectly perpendicular to the 
literature measurement for the disk angle PAdisk = 156°0 ± 0?3 
(Fig. Handel). 

Theoretically, an estimation of the inclination angle / stal is 
also possible from spatially and spectrally resolved interfer- 
ometric observ ations by using rapidly rotating stellar photo- 
spheric models (Domici ano de Souza et alj |2002). However, the 
limited accuracy and spectral resolution of our data does not 
allow the stellar inclination to be constrained in a meaning- 
ful way. We note that Fomalhaut has been observed intensively 
in broad-band interferometry with the yiNCI /VLTI instrument 
dDi Folco et al J 120041: iLe Bouquin et alj|2006l) . but these obser- 
vations, even when considered as a whole, did not allow to re- 
cover the inclination angle, mostly due to the very small ex- 
pected oblateness ratio of the Fomalhaut photosphere (~2%, 
O. Absil et al., in preparation). Weak constraints on the stel- 
lar inclination can nevertheless be extracted from simple phys- 
ical assumptions. The maximum permissible break-up veloc- 
ity for Fomalhaut is close to 390 km/s according to von Zeipel 
equations with typical values for A5V stars (2 M Q and 1.7Rq). 
Assuming a projected rotational velocity of 93 km/s dRover et al.l 
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Fig. 3. AMBER spectro-astrometric positions p(A) in the continuum (a) and across the Br-y absorption line (b). Colors refer to 
the wavelength bin, as shown in Fig. [2] The signature of the rotati ng photosphere (c ) is clearly detected and is compared to the 
debris disk and the planetary companion (d) imaged in the visible by Kal as et all (120081) . For the sake of clarity, the astrometric error 
ellipses are represented by their projection into the North and East directions. 



12001) . it gives a lower limit i stm > 15°. We conclude that the cur- 
rent existing constraints are compatible with the literature value 
for the disk inclination (/disk = 65.9°), but does not represent a 
significant test. Interferometric observations at higher precision, 
higher spectral resolution, and ideally with longer baselines may 
certainly answer this important point. 

Besides constraining the orientation of the stellar photo- 
sphere, spatially and spectrally resolved observations also con- 
strain the direction of its spin vector. Taking into account that 
the Br-y line is in absorption, we conclude that the South-East 
part of the star is moving towards us. If we assume that the or- 
bital angular momentum of the planetary companion points to- 
wards the same direction as the stellar spin, the Western side 
of the debris disk is located on the observer's side of the sky 
plane. Combined with th e observed increa sed brightness in the 
Eastern part of the disk (Kal as et ail l2005). it suggests that the 
dust grains are mostly backward-scattering. This finding is in 
contradiction with the well-known forward-scattering properties 
of circumstellar dust grains in our Sol ar system, generally as- 
sumed to be true in all debris disks (e.g., Weinberger et al. 1999; 
Kalas et al. 2005). Unfortunately, this question only appeared as 
a very interesting by-product at the time of the study, and our ob- 
serving strategy has not been specifically designed to answer it. 
Especially, no check star is available in the data-set to perfectly 
secure the sign of the AMBER phase. 

We did our best to calibrate the phase sign a-posteriori. 
First we checked that a positive delay corresponds to a nega- 
tive phase as measured by the AMBER reduction package we 
used. Secondly, we converted the phase curvature across the K- 
band (due to atmospheric refraction) into actual position in sky. 
We found that the blue part of the band is indeed shifted to- 
ward the zenith (as it should be). However, even if we are confi- 
dent, we cannot draw definite conclusions before a real spectro- 
astrometric reference has been observed. 



4. Conclusions 

From the technical point of view, these observations are 
a strong validation of the remarkable potential of the 
spectro-interferometer AMBER/VLTI, especially in its Medium 
Resolution mode and associated with the fringe-tracker FINITO. 
Stacking 17min of effective integration on a K = 0.94 mag star, 
we were able to reach a spectro-astrometric accuracy of +3 ^as. 



Associated with a spectral resolution/? = 1500, this was enough 
to resolve the rotating photosphere of Fomalhaut, with an appar- 
ent diameter of 2.2 mas and a rotational velocity of 93 km/s. We 
believe this to be the first detection of the astrometric displace- 
ment created by an absorption line in a fast rotating photosphere. 

Fitting our complete dataset, we find the position angle of 
Fomalhaut rotation axis to be PA star = 65° + 3°, perfectly 
perpendicular to the literature measurement for the disk angle 
PAdisk = 156°0 + 0°3. This is the first time such a test could be 
performed outside the Solar system for a non-eclipsing system. 
However, as when using the Rossiter-McLaughlin method, the 
inclination angle remains unknown. We can only conclude that 
there is strong evidence, but no definite proof, that the plane- 
tary system of Fomalhaut is in the equatorial plane of the central 
star. Additionally, by determining the direction of the stellar spin 
vector, we demonstrate that the technique is able to constrain the 
scattering properties of the dust grains surrounding Fomalhaut. 
Such measurements as presented in this paper are only possible 
thanks to the remarkable combination of a fast stellar rotation, a 
large stellar apparent diameter, the presence of a resolved disk, 
and of a planetary companion with resolved orbital motion. 

We plan to continue observing the Fomalhaut system to pro- 
vide even better constraints on the rotation of the central star, 
as well as to extend the study to other fast rotating stars hosting 
resolved debris disks. 
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